Examples 



We will give several examples of Green's formula. The first we will 
'build from scratch' so you get a sense of how this formula arises naturally. 
The last example shows how Green's formula works for a system driven at 
its resonant frequency. 

Example 1. The build up of a pollutant in a lake 

Every good formula deserves a particularly illuminating example, and per- 
haps the following will serve for the convolution integral. It is also illus- 
trated by the Mathlet Convolution : Accumulation. 

Problem: We have a lake, and a pollutant is being dumped into it, at a 
certain variable rate f{t). This pollutant degrades exponentially over time. 
If the lake begins at time zero with no pollutant, how much is in the lake at 

time t > 0? 

Solution. Let x{t)he the amount of pollutant in the lake at time t and a be 
the decay constant. For exponential decay we know that if a quantity p of 
pollutant is dropped in the lake at time Tjt then at a later time t it will have 
been reduced to the amount 

pg-a(T-Tt)_ (J) 

Here t — Xk is the time elapsed between when the pollutant is added and 
when we check how much of it is left. 

In our system pollutant is not being added all at once. Rather, it is drip- 
ping continuously into the lake. We break the interval [0, t] into n smaU 
pieces of width At as shown. 

Ar At At ^ ^ 

= To Tk Wi' Tn = t 



Let pk be the amount of pollutant added in the interval [t^, t^+i]. Since 
At is small we get the approximation 

Pk « /(Tfc)AT. 

(Remember /(t) is a rate; to get a quantity you must multiple by time.) 
According to equation (1) the amount of this left at time t is approximately 

«/(Tfc)ATe-''('-^^). 
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This is approximately the contribution to x{t) from the interval [r^, Tjt+i] . To 
determine the x{t) we simply siun up the contributions of all the intervals. 

x{t) « poe-<*-^''^ + pie-«(*-^i) + . . . + p„_ie-«(*-^«-i) 

~ (^f{To)e-<*-^°^ +/(Ti)e-<*-^i) + . . . +/(T„_i)e-<*-T"-i) j Ax. 

This is a Riemann siun. Taking the limit as At we get the convolution 
integral 

x{t) = f f{x)e-<^-^Ux. (2) 

Example 2. In example 1 we constructed our formula by slicing an interval 
into pieces. You should know how to do this. But, we prove theorems and 

formulas to avoid always going back to first principles. In this example we 
will solve the problem in example 1 using the differential for exponential 
decay and finding its weight function. (Of course, this DE was found by 
slicing an interval into pieces ) 

The DE with rest IC is 

k + ax = f{t), x(0")=0 

Its weight function w{t) is for f < 0, and for f > it is the solution to the 
IVP 

w + aw = Q, w{Q)=l. 
We get w{t) = e-«* for t > 0. 

Using Green's formula we again get the convolution integral (2) 
Example 3. Resonance 

Use Green's formula to solve the DE with rest inital conditions 

2x + 8x = cos{2t), x(0") = 0, x(0") = 

For t > 0, the weight function is the solution to 

2w + 8w = 0, w{0) =0,w{0) = 1/2. 
1 

The solution is w{t) = - sin(2f). For t > Green's formula gives 

ft 1 
x(t) = / sin(2T)-cos(2(f - T))dT. 
Jo 4 

2 
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This is an easy integral, we sketch the algebra to compute it. It uses the 

. ., . . / .N sin(A + B) + sin(A - B) 
trigonometric identity: sm(Ajcos(B) = — ^ ^ -. 

x{t) = /oSin(2T)icos(2(i-T))dT. 

= I J* (sin(20 + sin(4T - 2^) dr 

' • l^t.\ cos(4T-2t)" 

Tsin(2f) ^—^ '- 



_ fsin(2f) 
8 ■ 

This is the answer we expected from our earlier work with the exponential 
response formula. 
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